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Abstract 

A Hamiltonian eikonal model for multiple production in high 
energy hadron-hadron collisions is presented and worked out with 
the aim of providing a simple frame for various different observ- 
ables. 

An important role is played by unitarity which is built in 
by construction in the Hamiltonian formulation. The eikonal 
approximation allows both a very effective simplification of the 
dynamics, and facilitates the discussion on the relevance of pos- 
sible spatial inhomogeneities of the hadrons. The model is in- 
tended to describe only the hard interaction of the constituents, 
the structure of the incoming hadrons and the final hadronization 
processes are outside the scope of the present investigation. 



1. Introduction 

A Hamiltonian model for the description of the multiple production process 
in hadron-hadron collision is presented with the main aim of bringing together 
different observables within a unique frame. A particular attention is given to those 
features of the inelastic processes that can give informations on the proton structure, 
among these the relation between the hard-inelastic cross section and the inclusive 
production rates. 

The physical ingredients of the model are the following: 

• Although the underlying theory should be ideally QCD there is a sharp distinc- 
tion between soft dynamics, which provides the binding of the partons in the 
hadrons and also the final hadronization of the shaken-off partons, and hard 
dynamics that causes the parton scattering. 

• Hard collisions gives a finite transverse momentum to the partons, which re- 
mains however small with respect to the typical longitudinal momentum. Hard 
rescattering is included, but not hard branching of the partons. 

• Discrete quantum numbers like spin and colour are not taken into account. 
The observable quantities that can be computed are the hard inelastic cross 

section, the inclusive cross section for the production of back-to-back pairs of par- 
tons, the cross section for double pair production, the multiplicity distribution, the 
backward forward correlation between the produced partons. The effect on the 
observables of the hard dynamics, described by the Hamiltonian, will be distin- 
guishable from the effect of the hadronic structure, which is parametrized in an 
independent way. 

In the second chapter, after a short reminder of the eikonal formalism [1], whose 
use is suggested by the kinematical conditions, the general features of the model 
are made precise by defining the interaction and by choosing a definite partonic 
description of the hadron. In the third chapter the expression for the inelastic 
cross section is derived. In the fourth chapter a few observables related to the 
production process are calculated. In the fifth chapter the possibility of a non 
uniform distribution of the partonic matter in the hadron is considered and the 
effects of this hypothesis on the inelastic cross section and on the production rates 
are worked out. Since the whole treatment deals with the transverse variables 
an exploration on the possibility of taking, somehow, into account the longitudinal 
degrees of freedom and of some simple related consequences is presented in the sixth 
chapter where also an initial discussion of more general forms of parton distribution 
inside the hadron is sketched. 

One can find lot of previous treatment sharing important analogies with the 
treatment presented here, both in the eikonal formulations for the multiple produc- 
tion and in some purely probabilistic descriptions of the collision processes, which 
however, providing a form of unitarization of the transition probabilities, are in 
their final answer, analogous to the present formulation [2]. The main advantage in 
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the actual approach with the Hamiltonian formalism is that unitarity is explicitly 
implemented in all the different steps of the calculation and its role, with respect 
to the different observables considered, can be always traced back. 

A preliminary version of this work was presented at the XXVIII International 
Symposium on Multiparticle Dynamics - Delphi (1998)[3\. 



2. General features 

2.1 A short reminder of the eikonal approximation 

The description and justification of the eikonal formalism in high energy scat- 
tering has been presented in a lot of papers, so there is no point in re-deriving it. 
Only some features that are relevant for the next exposition are here briefly recalled 
[1,4]- 

The relative motion of two very fast colliding particles in their cm. frame is 
described by the free Hamiltonian; 

H Q = v • p + M/7 

One then adds an interaction term V which is for the moment left unspecified, 
but for the fact that it depends on the relative coordinate r = (B,z). At very 
high energy the speed remains practically constant even for sizable changes in the 
momentum, so the solution of this Hamiltonian problem is given by the wave- 
function: 

* (r) = (2^ exp h^ S- dz ' v ^ z>) + H (2J) 

The exponential factor yields, in the limit z — > oo, the S-matrix at fixed B. ^From 
this formula it appears clearly that the most important feature of the potential is its 
dependence on B whereas its longitudinal variables are always integrated; moreover 
the integration in dz 1 / (3 is equivalent to an integration over time. 

The considerations here presented are of pure kinematical origin, they keep 
their validity also when a more complicated structure is foreseen both for the inter- 
action term and for the incoming states as it will be done in the next section. The 
main improvement will be the description of the incoming states as systems with 
internal degrees of freedom [5,6] and thus the introduction of a set of transverse 
coordinates, one for each component. 

2.2 Description of the model 

The model describes the hadrons as sets of bound partons which, due to the 
interaction, may become finally free and eventually may be detected as jets in the 
final state; the hadronization process is not described. 
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The only detailed kinematics is the transverse one, the longitudinal is in some 
way integrated over, so also the longitudinal relative motion of the hadrons, which 
appears in eq.(2.1) as e ipz , is not explicitly written out. Since we deal with very-high 
energy collisions there is a sharp distinction between backward and forward degrees 
of freedom. We call a&, a/ the operators of the bound backward and forward partons 
and CbjCf the operators of the free partons. They have the standard commutation 
relations, every backward operator commutes with every forward operator and every 
a commutes with every c;they are local in b, the transverse impact parameter of the 
parton, this is possible since the size of the region relevant for the hard scattering 
is much smaller than the hadron size; so we can write a free Hamiltonian: 



The interaction that we want to describe is the hard collision of two bound partons 
that give rise to two free partons is such a way however that they keep their property 
of being either backward or forward. Thus the interaction Hamiltonian is written 

as: 



With this choice the interaction and the free Hamiltonian commute: [H Q , = , 
the theory however is not trivial, even though it has been simplified, the S-matrix 
can be written out in the form S = exp[— iHr] where r is an interaction time. 

2.3 Discretization 

The complete locality of the interaction in the transverse coordinates is both 
unrealistic and sometimes inconvenient, we consider an alternative with a finite size 
A of the hard interaction and a discretization of the transverse plane. The size A 
is related to the cut-off in the transverse momentum which must be put in order to 
be allowed to perform perturbative calculations, so the natural choice is A ps p^ 2 ; 
this choice leads also to the interpretation of r ps ~ \/A. The commutation 

relations become [Ay^A^J = Sij5 UjV , and so on. So in this discrete version the 
emission and absorption operators are dimensionless; also the coupling constant g 
is dimensionless, it is related to the previous coupling constant by A = gVK. The 
parameter O plays no role in the next treatment, it might also coincide with the 



H = " / d 2 b[al(b)a v (b) + 4(b)c„(b)] 



(2.2) 



v=f,b 




(2.3) 
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previously introduced u>. In this way we get 

v,j 

H I =(g/^)J2 H ^-H f , J (2-4) 

3 

H v j =CljA v j + AljCyj 
Also the S-matrix becomes discretized and it takes the form 

S = Y[Sj , Sj = exp[-i(^/v / A)r^, J • H u \ , (2.5) 

3 

With the previous interpretation of r, one obtains the simpler expression: 

S, = exp[-igH bJ • H fJ ] . (2.5') 

In order to apply this model one must choose a definite initial state; it will be 
factorized in the same way as the S-matrix: as far as its structure in a site j 
is concerned there are no strong indications. A possible choice, related to some 
theoretical ideas about the non perturbative partonic structure of the hadron[7,8], 
is the coherent state, so we may write 

|I>=IIl J> i > lI^eM-HlFbf + lFff^eMFbAl + FfA^l^ (2.6) 

3 

The vacuum state is also formally factorized; it has to be noted that the weight F 
of the coherent state may vary from site to site. For simplicity the index j will be 
not written out, whenever possible. 

It is useful to introduce the auxiliary operators through the definitions 



P = (C + A)/V2 Q = (C-A)/V2 
A=(P-Q)/V2 C=(P + Q)/V2 



(2.7) 



in this way we get 

= E { - l ,, r..i'\i + QtjQvj] , H VJ = P^P VJ - QtjQvJ (2-8) 

and it is also easy to express \I >j in term of the basis generated by P and Q. 



3. Inelastic cross section 
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3. 1 Contribution of a discrete site 

The first observable that can be calculated within the model is the inelastic 
cross section; it will be calculated using the discrete formulation for the states and 
for the S-matrix. In the basis generated by the operators P and Q the operator Sj 
is diagonal, so it is easy to calculate the matrix element Sj = j < I\Sj\I >j, it has 
the expression: 

S i= N2 E fc ,J fc , fc , (£1^ 

(3.1) 

The indices fci, &2, k$, k± refer to the quanta created respectively by Qb, Pf, Qf- 
The normalizing factor is given in eq.(2.6), actually 

iV = exp[-i(|F 6 | 2 + |F / | 2 )] . (3.1') 

By using the following representation, with a(3 = g, 

«p[- W *. - *,)(*, - M = (2.)-' / *«fc«p[i» + i Q u(fe - b) + i/M*3 - *.)] 

the multiple sum in the expression of Sj can be transformed into an integral. We 
use the positions T v = \F V \ 2 and we obtain: 

^3 = 2tt / dudv exp[iuv] exp[— T&(1 — cos aw) — T/(l — cos/?v)] . (3.2) 



Continuum limit 

When the distribution functions F v do not vary strongly from site to site one 
can devise a continuum limit * so that the partonic structure of the colliding hadrons 
is described as: 



/>=exp J rf 2 6[[-i[|/ b (b)| 2 + |/ / (b-B)| 2 ] + [/ b (b)a;(b) + / / (b-B)4(b)] 



> 



Here B denotes the relative impact parameter of the two hadrons. We consider the 
natural relation F m fy/A; therefore, if / is not singular, in the expression for Sj 
the factors T become small, so the exponential in the integral representation can 
be expanded and integrated term by term giving as a result: 

Sj?zl- T b T f (l - cos g) + \T h T f {T h + 2»(1 - cos g) 2 + ... (3.3) 



* This means precisely that the distributions are smooth, not that we consider 
the unrealistic limit p± — > 00 
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With our normalization the inelastic cross section at fixed hadronic impact param- 
eter B is 

a(B) = 2 < 7|(1 - US)\I > -|< 7|(1 - S)\I >\ 2 . (3.4) 

The product of the matrix elements S = ■ Sj is, as usual, calculated through the 
sum of the logarithms and the sum A Y^j is finally converted into the integration 
J d 2 b. The final result is: 

(j(B) = l-exp / d 2 b\-at b (b)t f (b-B) + 

J 1 n (3.5) 

\a 2 t b (h)t f (h - B) (4(b) +tf(b — B)) + • 



Two definition are introduced, a = 2A(1 — cosg), since this is precisely the param- 
eter which has the role of elementary partonic cross section, and t v (h) = |/„(b)| 2 , 
giving the transverse density of bound partons. 

The form of the inelastic cross section is quite usual. The second term in the 
argument of the exponential represents the rescattering corrections, which will be 
discussed below, where a nonuniform model of the hadron will be explored. The 
cross section arises from the integration over the impact parameter, so the result 
depends to a large extent on the properties of the density functions t v . 

As a simplest case the distribution can be taken to be completely uniform in b: 

t b (b) = p b 0{R-\b\) , t f (b) = p f #(R-\b\); (3.6) 

elementary geometrical considerations give \B\ = 2i?cos|7 with < 7 < w. The 
exponent in the integrand is given by the partial superposition of the two disks, the 
superposition area is 

W = i2 2 (7-siii7) . (3.7) 
and the cross section is expressed as: 
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J rf7sin 7 [l-exp[-^]] , (3.8) 



where v = ap b pfTiR 2 and £ = (7 — sin7)/7r = W/(7iR 2 ). 

With the previous interpretation of a, the numerical constant z/£ is the mean 
number of partonic interactions. In the limit of wholly uniform density the intro- 
duction of corrective terms like those appearing in the expression of Sj, eq. (3.3), 
amounts simply to a redefinition of v. 

It is possible to give a simple analytical form for Oi n in the two limiting situa- 
tions of very small or very large v. In the first case one gets 

(T in = kR 2 ■ v . (3.8') 
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In the second case we can start from the expression 

a in = 2tvR 2 [2 - D{y)\ (3.8") 

where the real function D{u) defined by eq(3.8) is monotonically decreasing, for 
large v it results D(y) m 2(6z/ 2 /7r 2 ) -1 / 3 ■ so that the geometrical limit of black 
disks AtxR 2 is approached. 

4. Inclusive cross sections and multiplicity distribution 

4-1 Pair and double-pair production 

The production of a pair is signaled both by the production of a backward 
parton and by the production of a forward parton, the Hamiltonian being fully 
symmetric. We make then the arbitrary choice of looking at the forward particles 
only; successively we shall investigate how much the rescattering processes may 
destroy the sharp correlation between backward and forward scattered partons. We 
start from the computation of the inclusive production from a single site * 

<X J >= J <I\S]c)C f S J \I> J (4.1) 

and we observe that can be treated as a number with respect to the forward 
operators, so the following relation holds: 

SjCfSj = C f cos(gH b ) + iA f sm{gH b ) . (4.2) 

The state \I >j is a coherent state of bound partons and contains no free parton so 
one obtains: 

< X 3 >= jKllA^Afsm^gH^I > 3 . (4.1') 

The matrix element in the previous expression is better calculated in the basis 
generated by the operators P and Q (details are given in the Appendix), and the 
result is 

< X 3 >= \T f [I - exp[-T 6 (l - cos(2<7)]] . (4.3) 

We can now go to the continuum limit, always under the hypothesis of smooth 
distributions t v (b); we find a problem in the presence of the function cos(2(7) instead 
of cos g, in fact the elementary partonic interaction enters in < Xj > through the 

* The production of a double pair in a single site is not included because the 
size A is defined by the hard interaction, so the production of four particle in a 
single site would not give rise, if not accidentally, to two distinct pairs, each with 
compensating momenta. 
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quantity k = |A[(1 — cos(2^)], which is related to a in this way: k = a - 1[1 + cos g\. 
The two constants coincide for small g, where both take the same value: a = k = 
g 2 A, higher powers in g 2 make however them different. Unitarity corrections at the 
level of parton-parton collisions are in fact different in the total and in the inclusive 
cross section. 

In the simple case of small values of T we expand the exponential of eq. (4.3) 
and we get the usual expression: 

X(B) = k J d 2 bt f (b)t b (b - B) (4.4) 

and when we consider the integration over the hadronic impact parameter 

D 1= [ X(B)d 2 B (4.5) 



the two integrals get factorized. When, however, in the same conditions we calculate 
the double pair production we do not end with a factorized expression (there are 
four factors and three integrations). In this case we obtain 

D 2 = k 2 J d 2 Bd 2 bd 2 b't f (b)t b (b - B)t f (b')t b (b' - B); . (4.6) 

In the limit of rigid disk, using the geometrical considerations and the definition of 
the previous section 3.2 the expression becomes: 

D 2 = 2(nR 2 ) 3 (p bPfK 2 ) 2 J sin 7 rf 7 = (nR 2 ) 3 (p b p fK ) 2 [l - 16/(3tt 2 )] . (4.7) 

A ratio of the quantities that have been now calculated and that is of phenomeno- 
logical interest is or e ff = [Di] 2 /D 2 , [9,10] which has the nice property, within this 
treatment, of being independent of k. It depends rather on the space behavior of 
t(b), namely on the hadron shape. Since in the rigid disk limit it results from eq 
(4.4,5) that 

D 1 = {nR 2 )K PfPb (4.5') 
the ratio we are looking for is given by: 

irR 2 

^ = — - w 2.2tvR 2 , (4.8) 

1-16/(3tt 2 ) ' v ; 

The previously calculated expression of Oi n is really the hard part of the total 
inelastic cross section, where as "hard" part we mean the contribution of all the 
events with at least one hard scattering. If we believe that, in going on with the 
total energy these events become dominating we would like to have this term not 
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too small with respect to the experimental <7j n , which in turn appears to be sizably 
larger, of about a factor 2, with respect to a e ff, so in this model we expect to 
approach at high energy the black- hadron limit which produces o~i n = AirR 2 . 



4-2 Multiplicity distribution 

The distribution of the multiplicities of the produced pairs is calculated by 
defining the projection operator over the number of free partons. Since there is a 
sharp distinction between backward and forward particles we can choose to take 
the forward parton as a signal of the pair production. For a fixed site j the number 
projector is * 

V n = - t : C^Ce~ c ' c : (4.9) 
n! 

The colon indicates the normal ordering of the C-operators, more precisely if the 
operators refer to the forward particles this is the forward normal ordering, and 
in no way it affects the backward operators. The properties of the V n operators 
are easily verified. They are evidently diagonal in the number basis and clearly 
< m\V n \m >= when m < n, for n < m = n+£ we get, through direct computation 

i-n i 1 u (ji + ey. (n + £)! M ^ . 

It could be more convenient to deal with the generator of the projectors: 



with Z =: e^ C : . (4.10) 



n\ V<9/i/ 

An auxiliary function is introduced 

Z(fj) =< I\S^ZS\I > (4.11) 

Calculations are strongly simplified by the normal ordering; through the same 
steps which lead from eq.(4.1) to eq.(4.3) one can get in fact 

Z(jj) =exp[i//T/] < I\exp[-lfjTfCOs(2gH b )]\I> (4.12) 

and one must remember that the functions of are not normal ordered. The matrix 
element of eq (4.12) is computed by expanding the operator into numerical Bessel 
functions and operatorial trigonometric functions, using the relation [11]: e zcosd = 



* Whenever possible the indices /, 6, j are suppressed 
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I (z) + 2 J2k Ik(z) cos(2k0), which allows to perform the same calculations leading 
from eq.(4.1) to eq.(4.3). The final expression is written as 



Z(p)=exp[^T f ] Io(-^T f ) 



+2 h{-\»T f ) exp[-T 6 (l - cos(2%))]] ' 



(4.12') 



According to eq (4.10) the derived multiplicity distribution, as seen from a single 
site is 



1 / d \ n 



n=-i 



When the production at the single site in not very strong it is natural to expand in 
the absorbing term A straightforward although a bit lengthy calculation yields, 
to the second order in TV 



K ™ =^ T /] n f 7 ^ 1 - W - cos(2^)]"exp(-|T / [l - cos(2s)] 
\T%[1 - cos(4^)]-exp(-iT / [l - cos(4</)]) 



+ 4- 



(4.13) 



The result is evidently a sum of two Poissonian distributions. The origin is the 
Poissonian distribution of the initial coherent state, which is modified in a well 
defined way by re-interactions. 

When re- inter actions are important there is no simple expression for the result, 
as usually we get a not too awkward expression in the extreme limit ,i.e., very large 
T&. In this case we could neglect in eq.(4.12') all the terms containing the negative 
exponent of T5. We translate [11] afterwards the expression containing the Bessel 
function into one expressed through the confluent hypergeometric function: 

and we obtain in this way for the multiplicity distribution the form: 



K n n -T f n exp[-T f ] ■ L n , L B = -(±) | M(±;n+ 1;2>) (4.14). 



The expression has been put into a form of a Poissonian distribution times another 
factor, this further factor factor L n is not a small correction however, it changes in 
an essential way the shape of the distribution. 



4-3 Forward-backward correlations 
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'X 



(4.16) 



Until now pair production has been described by looking at the production of 
a definite component of the pair. When the re-interactions are not very important 
this attitude is justified. If however there is a strong localized production in one 
site further investigations are needed. In this last case we can calculate within the 
model the variance and the covariance of the number of emitted partons. We start 
from the computation of the dispersion in the inclusive production in a single site 

3- 

< x 2 >=< I\S^C}C f C}CfS\I > . (4.15) 

The calculation goes along the same patterns as in the previous calculation of 
< X >. The result is 

< X 2 >= T] < I\ sm 4 (gH b )\I > +T f < I\ sin 2 (gH b )\I > 

and for the variance one obtains: 

E f = < X 2 > - < X > 2 = |T/[1 - exp[-T 6 (l - cos(2#)]] + |T|: 
(l - 2exp[-2T / [l - cos(2<7)]] +exp[-T / [l - cos(4^)]]) . 

The starting point to calculate the covariance is 

< W >=< I\S^C)C f ClC b S\I > , (4.17) 

Because of the coherent state structure of \I > one can show that 

S*C b S\I >= iF b sm(gH f )\I > . 

By using this relation the expression of < W > gets simplified to some extent. With 
some more work one recognizes also that 

Af sin(gHf)\I >=Ff cos gsm(gHf)\I > 
Cf sm(gHf)\I >=F f sin g cos(gHf)\I > 

From now on the calculations uses the results already seen, like eq .(4. 1,1', 3) and 
yields for the covariance £^5 =< W > — < Xf >< X b > the following expression: 

S /)6 = \T f T b sin 2 g [exp[-T & {l - cos{2g)}} + exp[-T/{l - cos(2^)}]] (4.18) 
In term of the quantities S f jb one can define the correlation coefficient as: 

w =p!^ < 419 > 
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A look to eq (4.16,18,19) shows, as expected, that for small values of T v , namely 
for small production, the correlation goes to 1, on the contrary when T v becomes 
large the correlation goes to zero. 



5. Non uniform hadrons 



5.1 Inelastic cross section 

We wish now to explore the possibility that the hadron, and its projection over 
the transverse plane, shows strong inhomogeneities in the matter density. This is 
represented by assuming the existence of black spots, that cover a limited amount 
of the transverse area, while a much fainter "gray" background smoothly fills the 
rest of the hadron. For the black spots the continuous limit (sect. 3.2) is not 
justified, but since they cover globally a small area, we can distinguish three possible 
kinds of collisions: The spot-spot collision, to be treated individually, the spot- 
background collision, where the background is treated as a continuum in b and 
finally the background-background which is nothing but the continuous limit already 
studied. In order to deal with the first case we start again from eq. (3.2), in the 
limit of very large Tj the integral will be expanded around the zeros of the exponent 
by setting: 

ua = 2nn + x > vj3 = 2Tvm + ; 
the subsequent Gaussian integrations over x an d 4> give: 



Sj « --^=^exp[i(27r) 2 mn/^] (5-1) 

9 V T f T b m,n 



In order to give an estimate of its value, the double sum is then converted into a 
double integration from — oo to +oo with the final result: 

Si " Sa = h^Tjr b (5 - 2) 

In the mixed spot-background collision, always taking eq (3.2) as starting point, 
it is possible to expand in one of the T v terms, keeping the full expression for the 
other one with the results: 

Sj!*Sf = l-T b (l-exp[-Tf(l-co8g)]) or 
Sj ^S b = 1 - Tf(l - exp[-T 6 (l - cos#)]) 

for the collision between a forward spot with the backward background and be- 
tween a backward spot and the forward background. The background-background 
contribution has been already given in eq.(3.3). 
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We are now in position to calculate the modification to the expression for 
the inelastic cross section, eq.(3.8), that are introduced by the hypothesis of the 
existence of the black spots. The gray distribution is again assumed to be uniform 
so the term in eq.(3.3) is constant and it will be denoted by S Q . The interaction area 
W is considered as composed by w small elements W = wA, the ratio £ = W/tiR 2 , 
will also be frequently used, so that w = ^(nR 2 /A). 

The quantum mechanical treatment of a nonuniform hadron is now presented 
in a form which contains certainly some rough simplifications that were unavoidable 
in order to deal with a system with an elevated degree of complexity. We may begin 
by assuming, at fixed j, a state of the kind: \I >j= \h >j \If >j 

\h>j= xexp[-i|F 6 | 2 ]exp[F b Aj] + yexp[-±|G 6 | 2 ] exp[G b Aj] |> (5.4) 

In this expression, and in the similar one for \If >, the terms F = l-Fle 1 ^ and G = 
\G\e lx denote two different thickness while the coefficients of the q.m. superposition 
are x, y. Since the two coherent states are not orthogonal the general form of the 
normalization condition is complicated: 

1 = \x\ 2 + \y\ 2 + 

exp[-i(|F| - |G|) 2 ] {xy* exp[|FG|[l - e^"*)]] + x*y exp[|FG|[l - e^*"^]] } 

When however the two thickness are very different the last term in the normalization 
condition is exponentially depressed and we are left with 

\x\ 2 + \y\ 2 1 . 

It is reasonable to expect that the same feature occurs also in calculating the rele- 
vant matrix elements, but in a significant case the calculations will be carried out 
explicitly so that the guess can be verified. In fact a non diagonal term of the matrix 
element of Sj is computed along the same lines yielding eq.(3.1,2) and the result is: 



Snd =W l j dudvexp[iuv] exp[-|(|F 6 | - \G b \) 2 - \{\F f \ - \G f \) 2 } 



x 



exp[-\F b G b \(l - e^- x)h cos au )] exp[-|F/G/|(l - e^~ x) f cos (3v)] 

(5.5) 

It is therefore enough to have one F very different from one G in order to obtain 
an exponentially suppressed contribution. 

One can see in this way that out of the 16 terms, that are produced in calcu- 
lating the complete matrix element of Sj, 12 i.e. the interference terms, are expo- 
nentially suppressed. 

It is useful to formalize this approximate result in the following way: let be 
\h >j= [|° > +H # >]b where 

< o|« >p»< «|o >p» and < o|o >ps< •]• >ps 1 (5.6) 
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(5.8) 



Then j < I b \I b >j~ 1 + A 4 2 an d the S-matrix element takes the form* 

Sj = i< °I+a* < •\]b,A< < > +H» >WI° > >]/,j( 1 +^ 2 )" 2 

which more explicitly gives 

Sj = [S + p 2 S b + p 2 S f + p 4 S s ](l + p 2 )- 2 . (5.7) 

In the limit of totally black spots, where we have the following values for the 
terms entering in eq.(5.7), see eq.s (3.3,3.6,5.3,5.2) 

S =l- ap b pfA 
S b =1 - p f A = S r] b 
S f =1 - p b A = SoVf 
S s =0 

the expression for the S-matrix is: 

s = n s i ■■ = s o ■ [i + v 2 v b + v 2 v f ] w (i + ^ 2 r 2w 

j . (5.9) 

=S Q S C 

where the factor S Q = S Q W gives the contribution of the background scattering. 

A much simpler expression can be written when fi is small and keeping into 
account that S Q is not very different from 1: 

S c = 1 - wp 2 A[p b (l - ap f ) + Pf (l- ap b )} (5.10) 

One may notice that the first factor, wp 2 A, represents the part covered by black 
spots within the interacting area of the hadrons at given B. A more general obser- 
vation is that the systematic use of the relations of eq. (5.5,6) eliminates the more 
complicated aspects induced by quantum mechanics and the answer is the same as 
if one would have taken a probabilistic distribution of spots in the transverse section 
of the interacting hadrons. 



5.2 Pair and double pair production 

The basic ingredient for these new calculations is always given by the expression 
in eq. (4.3), which must be particularized for the situations under examinations. 



* The parameter p measures the ratio of the area covered by the spots to the 
background; it is assumed small and real because the relative phases have no role 
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For simplicity and, more, in order to put in evidence the features of this particular 
model the background will be taken as thin so that in eq. (4.3) the first term of the 
expansion of the exponential is enough while in the case of the spots the exponential 
will be considered totally absorbing. The local production amplitude is the sum of 
4 terms, it will be indicated as 

<Xj >= [X , + p 2 X s , + p 2 X , s + n 4 X a , 8 ](l + p 2 )~ 2 (5.11) . 

The first term represents the pure background interaction, it is given by eq.(4.4), the 
second term represents a forward spot interacting with the backward background 
so Tf is large while Tj, = £&A, the third term represents the opposite situations, 
note that the resulting expression is not symmetrical because we are looking to the 
forward produced particles, finally the fourth term gives the effect of the spot-spot 
interaction. The four terms must be summed over the allowed values of j, i.e. over 
the position included in the interaction region. The sums may be expressed as : 

J^X 0i0 = K p f p b W 

j 

J2 X s,o =K Pb (T f /A)W 
3 (5.12) 

j 

J2x 8 ,s=UT f /A)W. 

j 

The dependence on the angular variable 7 is the same for the four terms, i.e. we 
can in any case write ^X = K^. The integration over the impact parameter is 
factorized, so we get for the inclusive production rates expressions as in eq.s (4. 5', 7). 
Actually 

D 1 = tvR 2 K , D 2 = -kR 2 K 2 [1 - 16/(3tt 2 )] 

The result for cr e ff is therefore precisely the same as in eq. (4.8). It has to be noticed 
that although Tf may be large (T//A) in finite and would remain finite even in a 
formal limit A — > 0. At first sight it seems that nothing is gained by introducing 
an inhomogeneity into the hadron, but in fact some news features are present. The 
expression of <r e ff is purely geometrical, it does not contains p, R is simply the radius 
of the area where the spots may be found; the expression of Oi n on the contrary 
contains dynamical parameters. It is therefore instructive to compare o\ n and o~ e ff. 

A clear, although unrealistic example is the limiting case in which the back- 
ground is so thin that it contributes negligibly to the inelastic cross section, then 
the S-matrix element, depends only on the spot-spot interaction and it has the form 

S = [l - ff 1 -ff 1" (5.13) 
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For small values of [i this expression may approach 1 and so the inelastic cross section 
becomes small, this situation would correspond to have few spots wholly black 
distributed in a wide and very thin background. The conclusion of this analysis is 
then that in order to have a e e < Oi n the hadron should be compact, i.e. without 
holes or transparent regions. There is however another possibility. The spots can 
in fact be distributed with some correlation among themselves, this possibility can 
be investigated, to some extent, within the model. The analysis is however formally 
awkward and not very conclusive, so it will be just mentioned but not reported 
explicitly. 



5. 3 Local effects of unitarity 



Until now the models of the hadron that have been considered in more detail 
cover two extreme situations, the case where there is not local rescattering and 
unitarity is relevant only to the whole hadronic interaction, allowing a description of 
the multiple disconnected partonic collisions, and the case where the re-interaction 
is so strong that, locally, the hadron is completely absorbing. We think useful 
to investigate briefly some intermediate situation, in this context sometimes the 
matter distribution in the hadron (the bound-parton distribution of the model) will 
be assumed to have a Gaussian shape, just in order to allow some explicit analytic 
calculation. 

We start from the inclusive productions: for the single inclusive we have in 
general in the model, comparing eq.s (4.3,4,5) 



D x = \ J t f (b + B) l-exp[-2/rf 6 (b)] 



d 2 bd 2 B 



(5.14) 



and for the double inclusive it results 



D 2 



= \ J tf(b 



+ B)t/(b' + B) l-exp[-2Kt b (b)] 



exp[ 



-2Kt b (b') 



d 2 bd 2 b'd 2 B 
(5.15) 

The densities t vary with the total energy of the process because, really, only the 
part of the partonic spectrum that can give rise to the hard scattering enters in 
eq.s(5.14,15) and this part certainly grows with the total energy of the collision. 
The simplest way in which this variation can be implemented is by rescaling the 
densities by a factor A growing with the energy. It is evident that in calculating 
cr e ff = [D 1 ] 2 /D 2 the factors affecting tf are eliminated in the ratio and the overall 
effect of A amounts formally to a rescaling of k. 

So, until 2nt(b) « 1 one expands the exponential and gets the expression 
already displayed which contains only geometrical elements, but when 2nt(b) m 1 
the final expression is less simple. It will be studied with the particular choice 
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£(b) = (h/tc) exp[— ^b 2 ] and with the definitions b = R + r/2 b' = R — r/2, so 
that 

Jt f (b + B)t/(b' + B)rf 2 S = (///2tt) exp[-i^r 2 ] = JF( r ) 

and 

1 - exp[-2Kt / 6 (R - r/2)]] d 2 lM 2 r 

When 2nt(b) < r and r is not too small, the integral can be separated into two 
parts, the first is approximately 

J F{r)[l -e- T ] 2 d 2 Rd 2 r 

with the bounds b < b , b' < b a , b 2 = ln(2 / uK/7rr)/ / u. The integration region in 
d 2 R has the same geometrical shape as the integration region yielding eq (3.7), only 
the role of the radius is played by b Q . The integration in d 2 r is always convergent and 
has no significant dependence on 6 G , so the conclusion is that this first addendum 
is proportional to b 2 and so to In k. The second addendum is certainly not growing 
with k, at fixed r, so at the end we get the result D2 oc b 2 oc Inn. The behavior 
of Di is much simpler to estimate, the double integral is factorized, the first factor 
does not depend on k, the second is proportional to hiK, for large k as it may be 
seen by explicit calculation. The conclusion of the analysis shows that the ratio 
giving cr e fj has a region where it stays essentially constant also if the hadron has 
no sharp boundaries but that at the end it will suffer a logarithmic increase whose 
origin has to be found in the local (i.e. at fixed b) unitarity. The moment at which 
these local unitarity effects begin to be relevant is when 2«t(b) is not too small with 
respect to 1. It must be remarked that the growth with energy is basically different 
for i7i n and for a e ^ since in the first case a logarithmic expansion is given already 
by the simple formula, which is "unitary" at hadronic level but not yet at partonic 
level, whereas in the second case the partonic unitarity is essential. The logarithmic 
growth is due to the particular choice of a Gaussian shape, an exponential shape is 
required to give rise to the square-logarithmic growth. 

6. On the longitudinal dynamics 

In the frame of the eikonal representation that has been used throughout the 
paper the longitudinal variables for the partons are in general expressed by means 
of the fraction of the longitudinal momentum carried by the partons. This is the 
natural choice for the bound partons, for a free parton the usual variable is the 
rapidity whose connection with the fractional momentum requires the introduction 



D 2 = \ 



F(t) l-exp[-2«t 6 (R + r/2)] 
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of the transverse motion. Since we are always interested in situations where the 
partons are tied in a clear way to the initial hadrons we shall use the fractional 
momentum x in every case, with the sharp distinction in forward and backward 
partons.* 

Starting from the formulation which is discrete in the transverse variables it is 
possible to introduce operators which depend on the fractional momentum and give 
the usual commutation relations like [A v j(x), A\ ^x')] = Si j j5 UjV 5(x — x') , and so 
on. and write the new forms of the Hamiltonians: 

Ho=J2[ dxU^Al ^A^ + Cl^C^x)] 

H I ={g/^)Y. H ^rH f , 3 . (6.1) 

3 

H v ,j= J dx[Cl 3 {x)A v ^x) + Al^C^^x)] 

The hadron momenta are denoted by II„ and the fractions x refer to the forward and 
backward total momenta according to the operator where they appear. In the same 
way the incoming states are described in terms of some distributions of longitudinal 
partons, the direct generalization of eq. (2.6), where the transverse factorization is 
maintained for the hadronic state \I >= Ylj \I >j and for the vacuum: 



| J >j= exp 



>3 



\ J dx(\F b (x)\ 2 + \F f (x)\ 2 )]+ J dx[F h (x)A\(x)+F f (x)A\(x)] 

(6.2) 

With this particular choice of the interaction term everything proceeds as before 
because it is again possible to define the auxiliary operators P(x) , Q(x). If we 
exclude black spots and go to the continuum limit the answer is: 

o-(B) = 1 - exp [-a J d 2 bdx t b (b, x)t f (b - B, x) + \a 2 ■ ■ ■] . (6.3) 

At this point, having put into the game the longitudinal variables it is clear that 
the choice t v = const finds no justification whatsoever; the most uniform choice 
corresponds to a uniform sphere, i.e. t v (h,x) = p v (x)\/R 2 — b 2 d(R — \b\) The role 
played before by the partial superposition of two disks, described by W — itR 2 £ is 
now played by the product of the two volumes that are superimposed 

U = J aAr 2 - bV-R 2 - (B - b) 2 $(R - \b\)ti(R -\B- b\)d 2 b . (6.4) 



* With this limitations it is not possible to introduce the hard production into 
the parton dynamics. It is possible to take into account re-interactions among the 
produced partons; in so doing the S-matrix element acquires a further phase. 
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This term can be written in a variety of alternative forms, but not completely by 
means of usual functions, a representation that turns out to be useful is: 



U = R 4 J e iX - B [j 1 (XR)/X] 2 d 2 X = 2nR 4 J ^(XB^XR^dX/ X . (6.4') 

The spherical Bessel function has the form ji(x) = x~ 2 sinx — x~ 1 cosx. The 
term U appears in the definition of the inelastic cross section and gives rise to a 
very cumbersome expression ,a much simpler form is found in the definition of the 
inclusive production and of the double inclusive production. The first term can be 
obtained also by direct integration: 

D 1 = (4tvR 3 /3) 2 k J p h {x)dx J p f (x)dx . (6.5) 

For the double pair production, the expression is: 

D 2 = (4ttR 3 ) 3 2R k J p b (x)dx J p f (x)dx * J [j^t^dt/t 3 . (6.6) 

The last integral gives: 

/•OO 

K = / [j^fdt/t 3 = 0.01433 

J o 

So we finally obtain: 



7. Conclusions 

The model presented and worked out in some detail allows a systematization 
of different aspects of the hard processes in multiparticle production and suggests 
also some interpretations in terms of hadron structure. The connection with QCD 
is not direct as it appears from the fact that the interaction term is quartic while 
the fundamental QCD interaction term is cubic. The reason is that in the actual 
approach the fundamental input is the parton hard collision, not the branching 
process. 

All unitarity corrections for the different processes are fully explicit, even 
though the final analytical expressions are given for the two extreme situations: 
weak perturbative corrections, very strong absorption. The analysis of the effects 
of inhomogeneities has been carried out purposely without specifying the possible 
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origin of this supposed property, in fact the quantities which have been considered 
depend on a global characteristic, the total amount of black area, and not on other 
details. The density distribution in the hadron is however relevant when a- lXi and 
cr e fi are compared. The representation of the hadron transverse area as a collection 
of elements of finite but much smaller extension requires a limitation of the mo- 
mentum transfer which cannot go below, say, of 5GeV so that the interaction area 
turn out of the order of (0.04 fm) 2 , small enough in comparison with the hadron 
extension. The assumption that the partons of the hadron build up, locally in b, a 
coherent state is a theoretical prejudice. The test is is difficult since the production 
process generally alters the multiplicity distribution, so that what finally one sees 
only a combination of hard and soft dynamics. To gain some better insight into the 
problem we have discussed, at the end of the 4th paragraph, the modification to 
the multiplicity distribution induced by hard rescatterings, 

Further questions can be raised: one is the possible coherence or correlation 
at different values of b,what could, at the end, also involve the role of the colour 
variables: in fact the transverse size of the hadron fraction over which the matrix 
element of S has been calculated is determined by the the transverse momentum 
transfer, which is a quantity relevant for the perturbative treatment. Larger sizes 
of transverse coherence may however show up in the multiple-production processes. 
Another question is why the incoming states should be locally coherent. More 
general initial states could be considered, but the new parameters one introduces 
in this way are at present out of control. 
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Appendix 

In this appendix we give some details on the calculation leading from eq (4.2) 
to eq (4.3) since this result is repeatedly used and it yields also a model for other 
similar calculations. Since the calculation is performed at fixed site, the index j 
will be omitted. ^From eq (4.1'), letting the operators Af act on the states \If >, 
it results furthermore. 

< X 3 >= T f < I b \ sm 2 (gH b )\I b > . (Al) 
Going to the basis generated by the P, Q operators it results also 
< I h \ sm 2 (gH b )\I b >=!-!(< I b \ exp[2ig(P^P b - Q\Q b )]\I h > +c.c.) . 
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The state is now expanded 

i J » >= 2( ^) 1 /2 m!n! i? r +n (n t ) m (-^) n i > 

m,n 

and the exponential acts now trivially on the Fock states so the final outcome is 
< I b \ si^( 9 H b )\I b >= \-\N*Y, 2 i m+ u)/2 mW \ F "r +n -M^irn - n)] + c.c. . 

m,n 

(A.2) 

A convenient way of computing the double summation of eq (A.2) is to perform 
first the finite sum at fixed I = m + n, and then the infinite sum over /. 

</ b |sin 2 (^ 6 )|/ 6 >= 

1 " iy2 E W |Fb|2 ' E (1) exp[2is(/ " n)} exp[-2i^n] 



i(l-exp[|F 6 2 |(l-cos(2^))]) . 

The last step requires the insertion of the actual value of the normalizing factor N, 
as it is given in eq (3.1'); then using eq (A.l) and the definition \F b 2 \ = T b we get 
eq (4.3). 
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